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•^^ . Abstract. The Jaynes-Cummings model, with and without the rotating wave approximation, 

is expressed in the conjugate variable representation and solved numerically by wave packet 
^/i , propagation. Both cases are then cast into systems of two coupled harmonic oscillators, 

reminiscent of coupled bound electronic potential curves of diatomic molecules. Using the 
knowledge of such models, this approach of the problem gives new insight of the dynamics. 
r^ . The effect of the rotating wave approximation is discussed. The collapse-revival phenomenon is 

O ,' especially analyzed in a non-standard manner. Extensions of the method is briefly mentioned 

in terms of a three-level atom and the Dicke model. 

1. Introduction 

The main concepts of a wave packet approach to cavity quantum electrodynamics (QED) where 
put forward in [i\. The idea is to formulate the model Hamiltonian, either the Jaynes-Cummings 
T^H- I (JC) model or the Rabi model [21 E], in terms of the quadrature operators for the field rather 

Tij" ■ than the commonly used boson ladder operators. The quadrature operators obey the regular 

^^ ! conjugate commutator relations like position and momentum of a particle. As such, the problem 

can be viewed as a wave packet evolving on two coupled and bound potential curves (arising 
1^^^ . due to the two-level structure of the considered atom). The state of the system is, of course, 

f— ^ I embedded in the wave packet, where roughly speaking the potential curves determines the 

internal state of the two-level atom and the vibrational states of the wave packet correspond to 
the field mode state. Thus, once the wave packet is obtained any quantity is easily calculated. 
K> ' The JC and Rabi models are related through the rotating wave approximation (RWA) in which 

j^ ■ rapidly oscillating terms are neglected in the Rabi model to give the JC one. Correspondingly, 

the JC model is exactly solvable and the physical quantities of interest are in general represented 
by infinite sums deriving from the quantized nature of the cavity field. These analytical results 
are valid only in the regime where the RWA is justified, and beyond these the Rabi model must be 
considered. In the conjugate variable representation, both models exhibit an avoided crossing 
between the two potential curves, crucially affecting their dynamics. Surprisingly, especially 
around the avoided crossing the solvable JC model shows a non-intuitive wave packet dynamics, 
while the Rabi model display a more expected evolution. 

In this paper we extend some of the results of pj and describe how the method can be 
generalized to multi-level systems. Explicitly, we deepen the analysis about the phenomenon 
of collapse-revivals and discuss the differences in the evolution between the two models. The 
approach to collapse-revivals in the JC and Rabi models is somewhat different from the one 
typically used for wave packets evolution in molecular and chemical physics. This is thoroughly 
considered and the relation between the two viewpoints is sorted out. In particular, what is 



characterized by the revival time Trev for the JC model is labeled the classical period in wave 
packet dynamics. Thus, the time scale set by Trev can become very long compared to what 
one would expect from the classical period, which comes about due to the internal two-level 
structure of the model. Unlike [l], this paper briefly considers multi-level systems, that is when 
the number of internal states exceeds 2. We first consider a simple example of a three-level 
A-atom coupled to a single cavity mode, which serves as a prototype of how to extend the 
theory. More interesting is the following system of A^ two-level atoms coupled to a quantized 
field, namely the Dicke model [4j. Here, the proper basis for the atomic subsystem is a collective 
one reminiscent of angular momentum states. Within this basis, the adiabatic diagonalization is 
straightforward and the 2^ potential curves are regained. An alternative approach is to use the 
Holstein-Primakoff representation [5] , in which the 2^ atomic (fermionic) degrees of freedom is 
replaced by a single bosonic degree of freedom. In such case, the coupled potential curves are 
represented by one potential surface in 2-D. 

We proceed as follows. In the next section [2] we introduce the model systems, the relation 
between the two models 12.11 is discussed and we give both models in the conjugate variable 
representation in 12.21 In this subsection we also mention how one may derive a semi-classical 
model that describes the population transfer between the two potential curves as the wave packet 
traverses the crossing. A deeper insight of the coupled dynamics is gained from the adiabatic 
representation, presented in 12. 3[ The following section [3] considers the non-intuitive evolution 
of the JC model compared to the Rabi model, while section U] in detail studies the collapse- 
revival phenomenon, emphasizing on parallels between these two models and other wave packet 
systems. Section [5] is devoted to generalizations of the method to multi-level systems. Firstly, 
in 15.11 the idea is sketched using a simple three-level system, and in 15.21 we consider the Dicke 
model. Finally, section [6] gives a summary and outlooks. 

2. The model system 

2.1. Relation between the Rabi and the Jaynes-Cummings models 

In most cavity QED experiments, their characteristics is well described by the Jaynes-Cummings 
model, in which single modes and atomic transitions are isolated and coherently coupled. This 
is achievable due to the strong atom-field coupling and the use of high-Q cavities and long- 
lived atomic states, such that losses can be discarded over typical interaction periods. Within 
the dipole approximation, a microscopical derivation leads to the Rabi model defined by the 
Hamiltonian 

H'Rabi = ^^ (a^a + ]^+ ^-(Jz + ng'o (a+ + a') (a^ + a) . (1) 
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Here, a^ (a) is the creation (annihilation) operator for the field mode; a\n) = \Jn -|- l|n -|- 1) 
(a|n) = ^Jn\n — 1)), the sigma operators are the standard Pauli matrices acting on the two-level 
atom; Uzli) = ±|ib), a"^|±) = |=f), Ox = <y^ ^ cr~ and Uy = —i{a^ — o"~), u> (Q) the field 
(atomic transition) frequency and q'q is the effective atom-field coupling. Before proceeding, we 
define characteristic time and energy scales by uj~^ and hco respectively, and hence introduce the 
dimensionless variables Hjiabi = H'j^^^^^/hoj, Q = Vl' /uj and go = ^q/w. The interaction in ([1]) is 
built up from four terms; a^a^ and a~ a corresponding to simultaneous excitation/deaxcitation of 
the atom and the field, and a^a and a^a~ originating from excitation of the atom by absorption 
of one photon and vice versa While the first terms are appropriately called non- energy conserving 
terms, the latter are termed energy conserving terms. In a rotating frame with respect to the 
first two terms of ([T]), the interaction constitutes precess with either the frequency |r2 -|- 1| or 
I O — 1 1 , and provided 

|0-1| < |17 + 1| (2) 



the rapidly oscillating terms may be left out. In this RWA limit, the Rabi Hamiltonian relaxes 
to the acclaimed Jaynes-Cummings model 



^ t i\ hn ^ ( + 

Hjc =(a'a+-l+ —0-2 + %o [cr^a + a 



(3) 



The applicability of the RWA is not solely determined by the requirement ([2]) , but also on the 
ratio go/^- For large values of go/^, the typical time scale of the interaction exceeds the one 
of free atomic evolution and the time-accumulated error arising from neglecting the non-energy 
conserving terms becomes important. Thus, apart from the condition ([2]) between the atomic 
transition and field frequencies, the validity of the RWA necessitates that go < ^ J6j. 

Given within the RWA, the number of excitations N = a)a + 2^z is a conserved quantity, 
and taking this symmetry into account the JC model is readily solvable with eigenstates 



\En)+ = COS ( - j |-F,n- 1) -Fsin ( - ) |-,n;. 



\En)- =sm\-\\+,n- 1) -cos 
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and A = 17 — 1. The corresponding eigenvalues read 



^n 
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together with the uncoupled ground state \Eq) = |— ,0) with Eq = —Q/2. 



2.2. Conjugate variable representation 

The algebraic approach applied to the models presented in the previous section is, in many cases, 
preferable to other methods. In terms of the Rabi Hamiltonian, no exact analytical solutions 
exist and approximate concepts are developed to find the required quantities in certain parameter 
regimes [T]. For example, this enables for perturbation expansions or truncation of continued 
fraction solutions. Non the less, the validity of the results are restricted to specific parameters, 
and to go beyond these one must consider numerical approaches. The procedure used here 
adopts the x-representation in which the boson operators define the conjugate variables 



P 



i—= [a^ 



7l("' + " 



(7) 



obeying [x^p] = i. In this representation, the Hamiltonians ([T]) and ([3|) become 
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It is convenient to unitarily rotate Huabi and Hjc by U 
transformed Hamiltonians 
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where Vh{x) = x^/2. In the current nomenclature both models are identified as two displaced 
coupled harmonic oscillators, in which, however, the amount of displacement and the character 
of the couplings are emphatically different. For x = the two potential curves possess an avoided 
crossing and the dynamics around such degenrate points have been thoroughly analyzed in the 
studies of excited diatomic molecules |7J. Typically, the crucial part of the evolution occurs 
when the wave packet transverses an avoided crossing. To a good approximation the potentials 
can be linearized in the vicinity of the crossing. Considering the wave packet as a classical 
point particle following the classical equations of motion, the population transfer between the 
two levels while passing the crossing in the Rabi model can be estimated by the Landau-Zener 
formula [H El [9] 

^/27rJ72\ 

8gov J ' 



Plz = 1 - exp 



(12) 



Here v is the classical velocity of the wave packet at the crossing, which, of course, depends on 
the initial state. Equation (I12p directly gives some indicatives of the system behaviours; for a 
large velocity v only a small fraction is transferred to the opposite state, and the same holds for 
a modest parameter VL which couples the two levels. In order for the wave packet to bisect the 
crossing region, its initial mean position Xi must satisfy Xj < or x? > 2gQ for an initial state 
starting out "on" the right shifted oscillator or Xj > or x? < —2gQ for the left oscillator. For 
such a state, starting out on a single potential curve, a classical estimate of the velocity at the 
crossing gives v ss J xj — (7q/2. The evolutions is said to be adiabatic if Plz ~ 1, diabatic in the 
opposite limit and mesobatic in the intermediate regime. The same arguments applied to the 
JC model is considerably less direct as the two potential curves are coupled by a "momentum" 
dependent term. This term is the background to the aberrance between the two models in a 
most unexpected way. For |J7| ^ goV^^ where n is the average number of photons (directly 
related to the velocity v), it follows that the equations can be decoupled, which corresponds 
to adiabatic elimination of the two atomic internal states. This will be confirmed in the next 
subsection. 



2.3. Adiabatic diagonalization 

The previous section introduced the idea of potential curves associated with the Rabi and JC 
models, and the concept of adiabaticity became somewhat clear in this picture. Here we will 
stress this even further in order to gain a deeper understanding of the dynamics. As the JC 
model is exactly solvable we only carry out the analysis for the Rabi model. 

The basis in which the Hamiltonian ([8]) is written might not be the most "optimal" one in 
the sense of decoupling the internal states. In most cases, the potential matrix describing the 
coupled dynamics has a smooth x-dependence, such that the characteristic length of derivative 
terms of the coupling elements decrease with the order of derivatives. A systematic way to 
remove the low-derivative terms from the off diagonals is by the adiabatic diagonalization. The 



unitary matrix Ui that diagonalizates a the potential matrix for the Rabi model is given in ^ 
by replacing the angle; tan (29) = 2\/2gQx/^. From the identity 



UlpUi = p — ayhdO, 



where 



ox 



(13) 
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it is clear that the transformed Hamiltonian is not diagonal. Explicitly we derive 
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The sizes of dO and d^9 measure the amount of adiabaticity [W\ lllj . From these we can draw 
several conclusions; a large favors an adiabatic evolution, so does a large x. The first is 
the regular adiabatic dispersive limit. The latter, however, is noticeably different from typical 
adiabaticicity in the JC model, but non the less it is intuitive since only close to the curve 
crossing is the adiabaticity assumed to break down. It seems that also $7 — > would give 
adiabatic evolution, which, however, is false since this limit describes the diabatic evolution. We 
recognize the adiabatic potentials 



X 



vTA^) = Y + ide{x)r ± \{x) 



(18) 



and the internal adiabatic basis states as | |) = cos{9)\+) + sin(0)|— ) and | |) = — sin(0)|+) + 
cos(0)|— ). The internal diabatic basis states are termed \u) = ^(|+) + | — )) and \d) = 

4=(|+) — |— )), with corresponding diabatic potentials V^ (x) = Vh{x it \/2go). Thus, a general 
adiabatic or diabatic state is is given by tp{x)\ | (|)) or ip{x)\u{d)) respectively. Examples of the 
adibataic (solid) and diabatic (dotted) potentials are presented in figure [TJ 

Labehng HJ^J = Had 
adiabatic diagonalization matrix U2 and so forth, leading to the expansion H^' 



UiHUl, the transformed Hamiltonian H^J defines a new 
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adiabatic Hamiltonian (jlSp shares great similarities to the JC one 
variables. 



We further note that the rather general form of the 

expressed in conjugate 



3. Regarding the p-dependent coupling 

The presence of the p-dependent term in the potential matrix of the JC model ([9]), makes its 
evolution, in some sense, more complex (or less intuitive) than for its companion the Rabi model. 
On the one hand, we know that the JC model is analytically solvable and any physical quantity 
can in principle be obtained from certain infinite sums. Non the less, the fact remains that the 
wave packet dynamics is still more involved due to this term. 
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Figure 1. This figure displays different examples of the adiabatic potentials V^^(a^) (solid lines) 
and diabatic potentials V^ (x) (dotted lines). 



The analysis will be restricted to either initial Fock number, or coherent states, which in 
x-representations read respectively 



V'n(a;,0) 
V'^(x,0) 



1 fl\^l^ 
, — H„(x)e 
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where H„ is the nth order Hermite polynomial and u is the amplitude of the coherent state; 
n = |i^p. Thus, in the case of a coherent state, the classical velocity v = J2n — 0^/2 at 
the crossing. For an initial Gaussian (J19p centered at Xi = 0, corresponding to field vacuum, 
one expects that the wave packet will split up and accelerate down towards the two potential 
minima [12]. This is indeed what we find for the Rabi model, while after all, this can not be 
the case for the JC model since we know that for an initial Fock state, the model exhibits Rabi 
oscillations. Thus, for the JC model the field wave packet will either Rabi oscillate between the 
two states iPq{x) and V'i(^) oi^ remain in iPq{x) throughout, whether the initial atomic state is 
1+) or |— ). This constrain of the wave packet to the origin is a direct effect of the p-dependent 
coupling, and holds for any Fock state ipnix) and therefore also for Fock states extending over 
large ranges x. The wave packet must be considered in its entirety as a single object, or in 
other words the coherence extending over the whole wave packet must be taken into account 
and one can not incoherently split up the wave packet in separate individual pieces. From this 
fact, understanding the influence of the p-dependent coupling becomes more challenging. None 
the less, we have that for a n photon Fock state Ap"^ = {p^) — (p)'^ = n + 1/2 and it is reasonable 
to expect that the "momentum" coupling will affect the dynamics due to the large spread in p. 



These conclusions are visualized in figure O which displays the squared absolute amplitude 
of the wave packet in the two models, together with the atomic inversion (a"z). We use the split 
operator method [13] to obtain the wave packet 
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and its squared absolute amplitude P{x,t) = |?/;+(x,t)p + \il)-{x,t)\'^ 
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Figure 2. The effect of the p-dependent coupling in the JC model. In the upper plot (a) we 
show the atomic inversion (o"^) between the two internal atomic states |ib) for the Rabi model 
(dashed) and the JC model (solid line) and the initial state is here; field in vacuum and the 
atom excited. The lower two plots presents the squared absolute amplitudes P(x, t) of the wave 
packet in the Rabi (b) and JC (c) models. The parameters are §0 = 2 and fi = 0.2. 



4. Collapse-revivals 

Collapses of physical variables are caused by dephasing between the constitute terms making 
up the quantum state. While the collapses are expected, more surprising is the phenomenon of 
revivals, occurring when the terms return back in phase. Hence, the evolution must be sufficiently 
coherent in order to be able to complentate the rephasing. It is clear that revivals in the models 
considered in this paper is a direct outcome of the quantized 'grainyness' of the field and thus 
a novel quantum effect [Ijj. In other fields, the existence of collapse-revivals has been greatly 
studied, and probably the most significant contribution owes the one of wave packet dynamics 
describing the vibrations in molecules [15]. Here the wave packet evolves in a bound electronic 
potential, where the discreteness derives from the vibrational eigen-modes in the particular 
electronic molecular state. For a harmonic potential, a wave packet bouncing back and forth 



in the potential, will reshape after one period of oscillation, which defines the classical period 
of motion. If the potential, however, is anharmonic the wave packet will not fully reshape after 
one classical period, bringing about the collapse. Depending on the degree of anharmonicity, 
the wave packet may reshape at later times characterizing the revivals. For a fairly localized 
excited wave packet, with average quantum number n ^ 1, we assume An/n <C 1 where An is 
the spread of quantum numbers. In this case we expand the eigenenergies accordingly, 

E{n) = E{no) + E' {n^){n - no) + ^^^(n - n^f + ^^^^{n - nof... , (21) 

where E'(no) = {dE{n)/dn)n=no aud so on. The various terms define different time scales 
according to 

rp ItV rp 27r rp 2?? /'OOA 

-^d — |_B'{no)|' -^r-ei; — |E//(„o)| ) -^ sup — \E"'(no)\^ V^^J 

characterizing classical, revival and superrevival times respectively. For a harmonic oscillator 
only the first of these term is non-zero and identifies the classical period Td = 27r/w. We note 
that, assuming zero detuning, A = 0, the JC energies expand as y^ = \/n + ^^{n — n) — 

Q-3/2 (^ — n)^ + ... . This is not, however, the standard way of deriving the revival times in the 
JC model. Normally one solves for the time it takes for consecutive Rabi frequencies to differ 
by a multiple of 2tt 

{2^n+i-2^n)Tl,, = 2Ti, (23) 



where Q.n = y — + don. The derived revival time is not identified with the previously defined 
Trev of equation ([22]) . but rather associates with the classical time T^. In other words, T^i does 
not correspond to the time for the constitute wave packets to bounce back and forth in the 
potential. However, in the JC model, as well as in the Rabi model, one has two internal degrees 
of freedom, and the dynamics is obtained from both the internal wave packets and their coupled 
motion. In particular, reshaping of the wave packet means that both internal wave packets must 
reform simultaneously. Thus, using the definitions (j22p . T^i is typically orders of magnitude 
larger in a multi-level system than for an internal structure-less wave packet evolution. Another 
way of picturing the phenomenon is that the constitute wave packets must overlap in phase space 
in order to give rise to interference, manifesting itself in the form of revivals. By introducing 
{x)± as the mean position of the two wave packets ip±{x,t), and similarly for the momentum 
{p)±, revivals occur when 6x = {x)+ — (x)_ =0 and 6p = (p)+ — (p)_ = simultaneously. In 
figure [3] we show the atomic inversion and in figure H] the quantities 6x and 6p as function of time 
t. In the x-direction, both models oscillates between ±10, while in the momentum direction the 
Rabi model has (p) exceeding 10 at some occasions. 

As argued, the revivals obtained in the Rabi and JC models correspond to the classical period 
Tci in equation ()22p . and not to T^ev arising from the anharmonities. Due to the internal two- 
level structure, the time Td can become rather long provided the proper adiabatic, diabatic 
or mesobatic potentials differ in their corresponding harmonic frequencies. However, if the 
frequencies almost coincide, Td does not need to be large as seen in the example of figures 
El (atomic inversion) and [6] (wave packet amplitudes). Here the parameters are such that the 
evolution is approximately diabatic and the two Rabi diabatic potential curves share almost 
identical curvature, causing the two wave packets to oscillate with nearly the same classical 
frequencies. The JC model shows more abstract dynamics, more reminiscent of the one seen in 
figures O and HI 

5. Extension to multi-level internal structure 

In m, the analysis was restricted to two- level atoms interacting with a single cavity mode. 
The generalization to multi-level atoms and/or multi-mode fields is straightforward. However, 




Figure 3. Inversion of the Rabi (black curve) and JC model (gray curve) for an initial coherent 
state with amplitude u = 7, and dimensionless parameters go = 0.15 and $7 = 5. 
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Figure 4. (Colour online) Inversion of the Rabi (red curve) and JC model (blue curve) for an 
initial coherent state with amplitude i^ = 7, and dimensionless parameters go = 0.15 and (7 = 5. 



the split operator method can only handle up to two dimensions, limiting the numerics to at 
most two modes, while the number of internal states may be as large as 20. However, other 
approximate wave packet methods exists where the dimension might be considerably larger [16j. 
Here we only consider extensions in the direction of multi-level internal states, and the many 
mode situation will be dealt with in the future projects. 



5.1. Three-level A- atom 

In the regular two-level atom model, losses of the excited level often affects the dynamics in an 
undesired manner. This can be circumvented by coupling two metastable "ground states" in a 
three-level A-atom configuration, adiabatically eliminating the excited state. However, the full 




Figure 5. Inversion of the Rabi (solid line) and JC model (dotted line) for an initial coherent 
state with amplitude u = 4, and dimensionless parameters go = 2 and Q = 0.2. 
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Figure 6. The squared amplitudes P{x, t) for the Rabi (a) and the JC (b) models corresponding 
to fig. El 



three-level dynamics show interesting features beyond the two-level atom situation [17j. The 
model studied here is used mostly to present the methods, as it can in fact be reduced into a 
two-level model. Additionally, generalizing this model to non reducible three-level systems is 
straightforward. 



For simplicity we assume the two lower atomic states, \gi) and \g2), to be degenerate, and 
further that they dipole couple to the excited state |e) through couplings Ai and A2. The 
Hamiltonian (without the RWA) reads 



Hlv = — + — + ^(Tee + [Al (o-gie + (^egx) + h (^-gae + (^egz) ]x, 

where aij = \i){j\- It may be written in the atomic bare basis {\gi), |e), I52)} as 
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Introducing the angles tan(6) = A1/A2 and tan(20) = 2Aox/0, with Aq 
potential matrix is diagonalized by 



U2 
As 8 is x-independent we find 

and one derives that 

The adiabatic potentials are 
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where the adiabatic state corresponding to Vq{x) is called the dark state which becomes 
degenerate with the state of V-{x) at x = 0. This degeneracy is lifted however, if a detuning is 
assumed between the two ground states. We note that the two adiabatic potentials V±{x) are 
similar to the ones of the Rabi model ()18p . Actually, the absence of a second non-zero diagonal 
element in the potential matrix gives the system a symmetry such that it can be separated into 
a 2 X 2 problem. The unitary transformation 
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casts the potential matrix in the form 



V\x) = U3V{x)ul 
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We may remark that like for the Rabi model, a semiclassical approach of this A-system 
naturally leads to the solvable generalized three-level Landau-Zener model [18]. Wave packet 
propagation methods in a related system has been considered in [19]. 



5.2. Dicke model 

Extending the Rabi Hamiltonian to contain N number of two-level atoms gives the Dicke model 

m 

2 2 

HD = ^ + ^ + nS, + -%S,x (32) 



with the total spin observables Sk = J2i=i ^k ^ k = x,z and a^ is the fc'th Pauli matrix for 



atom i. The scaling of the atom- field coupling by \/N is inserted to assure a well defined 
thermodynamic limit as A^ — > cxd. Operators acting on different atoms commutes, so that the 
adiabatic atomic ground state reads 

|i?o) = {cos(0)|+)-sin(0)|-)f^, (33) 



where 0N means direct product and as before tan(20) = 2gQx/\/NQ,. Formally, the adiabatic 
diagonalization procedure follows like for the single atom case, replacing Pauli matrices by total 
spin operators. We unitary transform the Hamiltonian by 



Ui = e-^'^'y (34) 



and the non- adiabatic corrections arise from equation (|T3l) substituting ay by Sy. From the 
theory of angular momentum we directly find the adiabatic potentials 




K^.(x) = y + m,y— + ^x2, m, = -N, -N + 1, ..., N - 1, N. (35) 

Here nis is the quantum number for the total spin in the z-direction. It is convenient to introduce 
the Dicke states \s,ms) being eigenstates of the total spin S\s,ms) = s{s + l)\s,7ns) and of the 
z-component Sz\s,ms) = ms\s,ms). The adiabatic ground state thus identifies with quantum 
numbers s = N and nis = —N, and applying the rotated annihilation operator J_|_ = f/l J+C/4 
repeated times to the ground state generates all the adiabatic states |iV, Tn^). The states with 
lower s quantum number can be obtained by standard methods. 

In the large A'^ limit, quantum fluctuations can in general be regarded as small and one may 
linearize these terms. The Holstein-Primakoff representation of the spin operators [5j has turned 
out to be an efficient approach in analyzing the Dicke model [20^ I21j. The spin operators are 
expressed in boson operators accordingly; 

S, = b^b.^, S+ = bWN -b^b, S- = Sl. (36) 

Before linearizing, we note from figure [T] that the low lying adiabatic potentials (we assume 
cold atoms and therefore a low temperature) either have one or two global minima, and this 
fact should be taken into account for when expanding in N~^. This is directly related to the 
presence of a quantum phase transition, with critical coupling g^ = vf^/2, between the normal 
phase of a vacuum field and all atoms deexcited and the superradient phase of a macroscopic 
field and atomic excitation [221 123] . Due to this it is convenient to coherently shift the boson 
operators as [20l [21] 

a^c + as, b^d + Ps, (37) 

where 



a.. 




go^N{l-f,^), /?, = ./ (1-^) (38) 



and 

i 1, 50 < g'i^ 

Thus, we note that in the case of a single global minimum the shifts vanish. Here c and d 
represents quantum fluctuations around the classical values ag and (ig- The resulting expanded 
Dicke Hamiltonian becomes 



For ^ = 1 (the ground adiabatic potential possess only a single minimum) the Hamiltonian 
becomes bi-linear and one may decouple the two boson modes into two disconnected harmonic 
oscillators. 

Using the Holstein-Primakoff representation we turn the fermionic degrees of freedom into a 
single bosonic degree of freedom. The cost of reducing the 2^^ potential curves to a single one is 
that the wave packet now evolves in 2-D rather than 1-D. Non the less, wave packet propagation 
in 2-D is easily performed with the split operator method, and will be consider in future works. 
Another algebraic method that could be considered is the "Schwinger's oscillator model of 
angular momentum" ^24], in which, however, two bosonic degrees of freedom is introduced 
instead of the fermionic subsystem. 

6. Conclusion 

The method of wave packet propagation, in which the cavity fields quadrature operators serve 
as conjugate variables, has been explored. It has been applied to the seminal JC model and 
its companion, the Rabi model, which is related to the JC one by the RWA. Various bases, 
different form the conventional bare and dressed bases, and their corresponding potential curves 
were consider. This new numerical approach is more commonly used in chemical and molecular 
physics, were it has been utilized for more than three decades [2S]. The effect of the RWA was 
discussed in this representation, and rather unexpected phenomena of the wave packet evolution 
appear once the RWA has been assumed. 

The main part of the analysis has been devoted to the collapse-revival effect of these models. 
Typically it is studied using algebraic methods, while we examine it from the shapes of the 
coupled potential curves. The terminology of collapse-revivals in wave packet models is not the 
same as the one for the JC model and we sort this out. The internal two- level structure of the 
Rabi and JC models may cause very long characteristic time scales, compared to the ones of a 
wave packet in a single anharmonic potential. 

Finally we have sketched how the method is extended to systems with more internal degrees 
of freedom, and here especially to the three-level A-atom and the Dicke model. For the Dicke 
model, we show how one may transform the internal degrees of freedom into a single " external" 
degree of freedom using the Holstein-Primakoff representation. More thorough research on 
multi-level and multi-mode systems are left for future works. 
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